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ABSTRACT 

This thesis examines the range error propagation and 
uncertainties associated with bistatic sonar operations. An 
equation for determining range to the target from the peeciyee | 
is explored for feasibility of practical applications. This 
particular equation does not require the separation distance 
between the source and receiver however, it does require an 
assumption of the mean sound speed over the two paths of the 
signal even though it could change drastically over a few 
miles. This thesis explores the contribution of eccentricity 
on the bistatic ellipse and the associated error. Examples 
demonstrate these effects by comparing cases of unequal mean 
sound speeds over the different paths at different values of 


eccentricity. 
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I. INTRODUCTION 
A. BACKGROUND 


Early ASW relied solely on active sonar, or Echo Ranging. 
A pulse of sound was projected into the ocean and if a 
submarine was in the path of the signal, an echo would return. 
The round-trip time of travel of this signal was used with the 
sound speed in the ocean to determine the range to the target. 
This method of ASW was very effective until submarines began 
to exploit the water conditions of the ‘Geenn and became 
evasive, i.e. , hiding in shadow zones, fronts and eddies, etc. 
The projected pulse from the active platform usually could be. 
detected by the submarine at a substantially greater range 
than that at which the active platform could detect the 
reflected pulse from the target because of the greater 
transmission loss (TL) over the round trip propagation path 
(Ref. 1]. This was one of the factors leading to interest in 
passive sonar. Additicnally, passive sonar: often utilized | 
lower frequencies having smaller losses, which could allow 
longer detection anges (Ref. 2). | | | | 

The structure of the speed-of-sound profile leads to a 
deep sound channel which can be exploited in convergence zone 
(CZ) propagation. 7 This propagation mode with its 
intrinsically lower acoustic losses provides longer standoff 
ranges, keeping the detector out of torpedo ranges of the 


submarine. On the other hand, a disadvantage of passive sonar 


oe 

















is that it requires a relatively time-consuming process of 


bearing rate analysis known as Tayget Motion Analysis (TMA) 


to determine range to a target [Ref. 3]. - Because of this, | 


active sonar operations are still important, particularly for 
iast minute, close-in prosecutions where exact range must be 
determined quickly for weapon deployment. | 


As technology advanced, passive sonar became more 


sensitive, ‘permitting longer detection ranges. Submarines, 


however, reduced their radiated noise levels nullifying the 
advantages of the more sensitive sensors. This progressive 
technology also made way for the continued improvement of 
_active sonar, extending capabilities to much lower 
frequencies, thus permitting active sonar exploitation of the 
deep sound channel. Bistatic Active Sonar (BAS) has brought 
about longer range aatsction by active sonar, and with it the 


advantage of immediate determination of the range to the 





target. Consequently, the vulnerability of active sonar 
platforms has been reduced with the’ technique or bhistatic 


sonar operations. 


-B. BISTATIC SONAR 


Conventional active sonar is monostatic, meaning 


hi 
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Pigure 1 Monostatic Arrangenent 


gonar sensors provide an estimate of the bearing of the target 
and the range to the target is determined indirectly by the 
equation | 


Range = = t-v (1.1) 


where V is the sound speed in the ocean, and t is the round 


trip travel time of the pulse and echo. | 


Bistatic sonar has a separated source and receiver on two 


different platforms. Figure 2 is a representative bistatic. 


arrangement [Ref. 4]. The receiver is positioned away from 
the source to receive the reflected signal from tke target. 
“In the bistatic case, the source, receiver, and target form 
a triangle which complicates the problem compared to the 
monostatic case because the range (R) is more difficult to 


determine. The source is presumably detached and away from 


the threat, reducing its own vulnerability. An active 




















Figure 2 Bistatic Arrangemet 


signal, or "ping", originates at the source, travels to the 
target, is reflected from the target, and then travels to the 
receiver as shown by the arrows. ‘The *~tal travel distance 
of this signal is X. The receiver also receives a direct 
signal from the source over the path S. The angle (y) is the 
separation angle of the source and the target ne:.sured at the 
receiver. ‘The ultimate goal of such an operation is to obtain 
the bearing and range of a target with respect to the 
Pacelver: Since the bearing is directly obtainable from the 
directional sonar receiver, determining the ands (R) is the 


problem that is to be treated here. 
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Ii. STATEMENT OF PROBLEM 





A. RALGE EQUATION | | =f 

An equation must be derived that solves for range (R) to - 
the target in terms of known quantities. For our purpcses, we 
shall assume that the total travel distances of the acoustic 


signal over the two paths between source and receiver, the 


distances X and S in Figure 2, are known. Then, given the 2 
bearing y, we have 7 
and, from the law of cosines, i 
go? = R2+S52-2RScosy. (2.2) _ :, 
Substitution of the first equation into the second alluws the . 
elimination of Co. ; 





(X-R)? = R2+57-2RScosy | (2.3) 
Expanding the term on the left gives 
X?-2RX.+ KR? = R2+S?-2RScosy. (2.4) 
Rearranging, we have | | | 
| 2RX ~.2RScosy a X2-97, . (2.5) 
Solving for R results in [Ref. 5] | | 


- __xX*-s? : | 
R " 2x-Scosy) * _. = 


‘The above equation gives R in terns of S, X, and Ys 
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The time of travel of a siqnal from the source directly to 
the receiver over the path S is t,, and the time of travel for 
the indirect route X is t,. The difference between t, and t, 
is 

t = ty-ty. | (2.7) 
These time definitions, combined with the mean sound speed. 


over the propagation: path, determine the distance that the 


signal has traveled. Use of 

x = Vt, and s = Vt, | (2.8) 
yields several different equations involving different 
combinations of the variables; V, S, X, ty, t,, t and y (Ref. 


6). One such combination can be derived as follows. 


Substitution of (2.9) into (2.6), results in 


(Vt,)?-(Ve,)? | 
ae 2 (Vt,-Vt,z0sy) ' eee 
and, the substitution of (2.7) for t, provides 
See eee 
R = (V(t+t,) )*-(ve,) (2.10) 


2 (V(t+t,) -Vt, cosy) ° 


Expanding both the numerator and the denominator leaves 


Ver+2veee + vie? - ye? 


R= ——————__*___*_..__{, 
2(Vt+Vt,-Vt, cosy) 


(2.11) 


and simplification gives 
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Vrt+2vetc, 
2V(t+t,-t,cosy).’ 


k= 


Now, use of 
1 - cosy = 2sin?t 
produces 


— Vt(Vt+2Vv¢t,) 


and therefore 


t(Vt+2Vve,) 
2(t+2t,sin?y/2)— 


Now, use of the relationship 


ve 2 
Cs 
gives 
| t St .5 Sts 
7 t, e. 


2(t¢+2t,sin?y/2) 
and finally 


St(.t+2¢,) 
2t,(t+2¢,sin?y/2)— 


2v(e+t,-t,+2¢,sin? (T)) 





(2.16) 





(2.17) 


(2.18) 


This equation, commonly encountered, is a function of S, 


t, t,, and y. Measured quantities of the variables used above | 


are often known to much more precision than the value of the | 
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mean sound speed V over many miles. Thus V has the largest 
amount of uncertainty, because the sound speed profile at the 
receiver must be assumed over the whole path of the signal. 
(Most bistatic sonar operations use a formula similar to 
Equation 2.18 to determine range from the receiver to the 
_ target (Ref. 6].) However, this thesis explores the use of. an 
| equation that contains V explicitly rather than S and also 
allows the treatment of different mean sound speeds. By using 
. the variable S, there is an inherent assumption that S has a 
readily obtainable value. To see the advantages using V 
rather than S, it is appropriate to consider the role of ASW 
Command Control and Communications (Cc ), few operational facts 


must be discussed. 


B. ? CONSIDERATIONS 

Since S ie the distance between the source and receiver, 
an accurate position of the source must be known for the 
calculation of the separation distance, S. This requires that 


reliable communications be maintained between the two vessels 


carrying source or receiver (Figure 3). snips that are within © 


"line of sight" usually communicate by UHF (Ultra High 
Frequency) radio, which is extremely reliable, but very 
limited in range. | When ships are separated by a distance 
‘qreates than “line of sight", approximately 20 nm, other means 
of communication must be utilized. One such method is the HF 


(High Frequency) radio which has a longer range because of the 


‘radio wave propagation that follows the curvature of the earth | 
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Figure 3 Communications Illustration 


(Ref. 7]. The dvgadvantede of HF sommini cee ions is that this 
type of communication is extremely dependent on atmospheric 
| conditions, since atmospheric "ducting" is required. Another 
disadvantage of HF communications is the vulnerability of the 
HF transmitter to direction finding (DF) (Ref 7 8) by an 
opponent. esp iaveipiles tharecare Soviet HFDF sites around the 
world that can triangulate HF, transmissions to determine the 
location of the transmitter. A third method of communicating 
over great distances is UHF Satellite communications (SATCOM) : 
In this method, the UHF transmission is "bounced" off of a 


satellite to another station over the horizon. SATCOM is 


extremely reliable, and is indeed a preferred method, but in | 














the ee of "killer satellites," dependency on satellites can 
be extremely limiting in operations. 

If (2.18) were being used, then a communication link must 
be maintained. If this link is broken and communications ars 
lost, there are two options to continue bistatic operations. 
The first option is that ‘the source and receiver move to 
within line-of-sight UHF range, typically 20 nm. This would 
reduce the advantage of the bistatic concept because the 
receiver is not "out front" exploiting the subma>ines 
reflected energy, and decreasing’ the source's vulnerability. 
The second option in case of lost communications is for the 
receiver ship to approximate the source's position by applying: 
the last known course and speed to the last known position, 
also known as Dead Reckoning (DR). This . introduces 
significant error in the variable Ss. 

An alternate means of locating the target using bistatic 
operations would be based on an equation that removes the 


necessity to know the distance S directly. 


Cc. AN ALTERNATIVE EQUATION 
To avoid the difficulties mentioned: in the previous 


section, we return to (2.6), 


oe X*-s? 
2 (X-Scosy) 


and use the relationships 


X= Ut, and $= Vt, (2.19) 
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to introduce differences in the mean sound speed over the 
aifferent paths, where U is the mean sound speed over the 
reflected path and V is the mean sound speed over the direct 
path. Subsiituting (2.19) into (2.6) results in 


re y2 
(Ut,)° - (VE,) (2.20) 


~ 2 (Ut, - Vt,Cosy) * 


and by dividing the numerator and’ denominator by Vv’, 9 


becomes 


vfs) ~ (te)? 
- 3 , (2.21) 
== tx ~ t,cosy| 


A scaled time, 1f,, can be introduced to deal with these 


differences in speed, where 
a a : 
Eugation (2.21) can now be rewritten as 


_ v(<2- ¢? 


; (2.23) 
. 2(t, — €,cosy) 


The equation still requires values for'r,, ts, and Y, but 

S is no longer directly used. This thesis examines. the 
| utility of this equation (2.23) as an alternate to (2.18). 

An important point to note is that in order for bistatic 

operations to take place, convergence zone conditions. must 

exist, and for these to exist there must a minimum of 200 


fathoms of depth excess for reliable propagation (Ref. 3}. In 
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essence, bistatic operations usually must take place in the 
deep ocean and away from typical frents and eddies where water 
conditions change drastically over short distances. 

The feasibility of using this equation will be examined by 
means of error analysis. We include in this analysis the 
variation in mean sound speed between the two different paths 
as well as the uncertainties in bearing and in times of flight 


or signal and echo. 
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III. DEVELOPMENT 07° ANALYSIS MODEL 


Analysis of (2.23), 


V(t2- 2) 


R* 3 (t,- t,cosy) ° 


is conducted from the point of view that it is the equation of 


an ellipse. Indeed, any bistatic sonar geometry can be 


described in terms of ellipses, the properties of which relate 


directly to the sonar parameters. 


A. INTRODUCTION OF ELLIPSE 

Figure 4 is a diagram of a bistatic arrangement with a 
superimposed ellipse, where a is half the length of the major 
xici’s (the semi-major axis), and S is the distance between the 


two foci. The eccentricity (e) of an ellipse can be defined by 


the ratio e=S/2a. Use of the previously-defined relationships 


and properties of ellipses then provide additional expressions 


Fog Sei es 


gee 2 oS. 
 X GtR Ut, ty; 


(3.2) 
‘Equation (3.1) reveals that the effect of a difference between 
the mean sound speeds U and -V over the two paths can be 
completely accounted for by making an aaj uatwant in the 
eccentricity e of the ellipse. The manipulations that follow 
will result in convenient expressions: for the fractional 


errors in range arising from uncertainties in timing errors, 


13. 
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differences between U and V, and bearing errors. The value of 
eccentricity may range between 0 and 1, where e = 0 is a 
perfect circle ande «iisa very flat ellipse approaching a 
straight line joining the foci. These elliptical 
relationships will be useful in describing the error 


associated with different geometries. 





Figure 4 Elliptical Arrangement 
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B. DIFFERENTIATING TERMS 
The error in range, AR, can be expressed in terms of the 


errors in V, 1,, t,, and y as 





AR = Rays OR 


av. Oty 


SR At + FE hts + ay, (3.2) 


or, in the form of a fractional error, as 


AR. 1 AR, 1 OP 1 1 OR 
R *R2ov"? RO, iS Roe, Be, 4 te* "ea 


(3.3) 


The coefficients of the error increments AV, Ar,, At,, and Ay 
on the right side of the equation may be obtained by 
logarithmic partial differentiation of (2.23). The 





coefficient of AV in the first term on the right is 

ae . 
“Ae? y’ (3.4) 
The coefficient in the second term is 


OR 2ty 1 


= - —__—___. : (3.5) 


1 
RO, ti-t2 t,-Ccosy 


The coefficient in the third term is . 








1 0R | ~2e, _ cosy : 
Rot, +t2-t2 ,-t,cosy- | (3.6) 


' Pinally the coefficient in.the y term is 
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(3.7) 


If these coefficients are assembled, (3.3) may be rewritten as 


AR ay | ated Jax 


Ro Vo |e2-c2 t,-tcosy 
a (3.8) 
2c Y  tesin 
ae Ss a cos At,;- Nf Y A 
ty-t,2 ty~ C cosy ty-t, Cosy 
Oy At At At At 
Note: t,=2t,, At,= ZAt,, —% = —%, and 8 = VA5s 
V V Ty ty ty U ct, 


By incorporating the relationships described in (3.1) the - 


above equation may be rewritten as 





AR, AV {2 .__ 1) Ate 
R V 1-e*. l-ecosy) ty, 


(3.9) 








+S -2e cosy |e _ _esiny Ay. 


e? ” Tecosy ty 1-ecosy 


C. MEAN SQUARE ERROR 

‘The fractional error in range can be obtained through the 
use of this equation if all errors in the variables V, r, ty, 
and y¥ are known completely ; both in magnitude and in 
direction. On the other hand, if ali errors. (i.e., 
uncertainties) in the variables Vv, 1,, t,, and y are assumed 


to be statistically independent of each other, they may be 
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combined into a sean square error (MSE). The MSE is used to 
determine the total combined error resulting from errors in 
variables in the range equation. It is calculated by taking 
the square root of the sum of the squares of the individual. 
errors. Each of the terms has an ‘uncertainty value (A) 
associated with that particular variable, i.e., At,. The 
correct numerical uncertainty values of each variable used 
operationally has a security classification and will not be 
used in this thesis. To simplify the notation in the model 
the coefficients of the error terms in (3.8) will ‘be expressed 


as Ay, As or A;, where 


Ae ® 2 ee | ° 
x \i-e2 1-ecosy)’: : | 
-2e, cosy ie : 7 ; 
1-e? l-ecosy) . 


5 ee siny 


1-ecosy 


A; = 


The use of these coefficients gives the fractional error in 


range in the form 


AR. AV.,: At At, | | 
‘AR = AV, 7 +A; S+A,AyY. (3.10) 
Ty . Ty ' : i a 


R V 








Ay 


The mean square fractional error is thus’ 


ct 














; 2 2 
5 (AR) wea (Ap) rae[ Se) as] +a? (Ay)?. (3-22) 


ir which it is assumed that the errors on the right side of 
the equation are statistically independent. Equation Sead 
calculates error in terms of the fractional timing errors 
referenced to a time base of t,. If it is preferable to 
express timing errors over each path referenced to the travel 


time over each path, then 





+A,AY. © (3.122) 
Cs 














The mean square fractional error is then changed to © 


AV Ar,\? At,\* | ss | 
Raa on 





* (AR) -2,* 


To understand the use of these equations ouErose we 


consider the general method by which the range from “=e 
rece.ver ship to the target (R) and the associated range error 
can bn found. Assume, for example, that the signal travel 
times t, and t,, and the azimuth angle y all have been 
measured. Assume ’ further, that the average sound speed V 


over the direct path from source to receiver is known. Also 
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let .U, the ‘average sound speed over the secondary path, 
s-urce-to-target-to-receiver, be the same as V, so that T,=t,. 
The range R can then be calculated from 2.23, where 


v(t? - ¢2) 


Assume also that the uncertainties in cime, At, and At,, 
and the uncertainties in azimuths, Ay, are all known. The 
fractional error (uncertainty) in range associated with each 
of the uncertainties in the separate variables can be obtained 
as follows. | | 

For ‘timing error 4a; t, alone the corresponding fractional 


error in range is 





R) = Ay P | (3.14) 


We obtain A, from the graphs of Fig. 5 through 13 after first 
calculating the eccentricity, | 


for this particular case being studied. ‘Similarly for error 
in t, alone use | 


’ 
’ ‘ 


At, 








[= }. Ay" = eA, (3-25) 


and for error in y use 
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(=4) = A,AY, (3.18) 


where Ay must be in units of radians. If more than one of the 
independent .variables is uncertain, then either (3.11) or 


(3.13) for the combined errors shou2d be used. 


D. ERROR ANALYSIS 

A computer program was written to study the behavior in 
range error in each of the three variables t,, t,, and y. The 
program, written in FORTRAN, was used for nine different 
values of eccentricity e of the basic range ellipse, from 0.1 
to 0.9. The output data of the program, the values of Ay, Ag, 
and A,, were plotted in the form of curves from which error 


estimates can be made for any reasonable scenarios. 
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IV. RESULTS AND VERIFICATION 


A. RESULTS OF MODEL 

Figures 5 through 13 are plots of the output data of the 
previously mentioned Fortran program. All of these plots 
indicate some relationships which will be examined further. 
Each of the nine figures is a graph of the error coefficients 
for a different eccentricity e, beginning with Figure 5 for 
the case of eccentricity 0.1 and continuing to Figure 13 with 
eccentricity of 0.9. In all of the plots the A, values are 
represented by solid lines, the A, values are represented by 
dotted lines and the A, values are indicated by dashed lines. 
The error coefficients (ordinates) give the magnitudes of the 
errors associated with each of the three terms. The baseline 
number (S) used for input is inconsequential because the value 
of e characterizes the driver in ‘the error propagation. Each 


of the three terms will be evaluated individually. 





The A, term is the indicator of the behavior of fractional _ ! | 
error associated only with uncertainty in the value of t. 
The graphs in the figures show how Ay and the other error 
coefficients vary with bearing angle Y- The graphs reveal 
that for low values of eccentricity e, the error coefficient oa a 
A, begins small and shows a small increase as y increases to 
180°. For large values of e, the values of A, remain small 


for low|values of y, but as y increases, the A, value | 
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increases rapidly until about mid ranae of y, approximately 
70°, when it levels off and apnroaches a maximum at 180°. 
This information indicates that for small values of e, the 
error contribution by the t, term is small for all values of 
y- On the other hand, for large values of eccentricity, the 
error contribution of t, begins small for angles less than > 
10°, but rapidly increases for larger angles until it steadies 
out. In practical terms, recall that a small value of e means 
a small separation distance (S )between the source and 
receiver compared to the distance between the receiver and the 
target (R). Consequently, a smaller S means a smaller error 
associated with t... For larger values of e, the error 
contributed by t, increases Crastically for separation angles 
(y) greater than 10°. 
The A, coefficient characterizes the effect of - 
uncertainty in the time .t,. Thesé values are represented by 
the dotted line on each of the nine graphs. For y of 0° and. 
| 180° this error coefficient is about the same as for the 
| previously determined A,. For small values of e, A gradually 
decreases to zero, at about 70°, and then increases until it 
has the Same value as A, at 180°. ‘For higher values of e, the 
| value of A, drops meee rapidly to zero at angles less than 10° 
and then approaches the values of A, as it increases sharply 
and behaves asymptotically done about 70°. The A, coefficient 
goes to zero at moderate values of y, about 70°, for low e, 


and it goes ‘to zero at small angles for large e. When the 
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source and receiver are closer to each other than the receiver 


ae to the target, this angle at which the error. vanishes is 
around 70° to 80°, and for large values of e, this angle is . 
much lower, near 10°. 

The error coefficient A,, associated with the y term, is 
zero for y values of 0° and 180°, as expected. For smaller 
values of e, the error coefficient increases slowly until it 
peaks at about 90°. For small values of e, the contribution 
of the y remains small compared to the other two terns. por: 
larger values of e,. the maximum value of A, increases 
significantly, arid it "shifts" to smaller angles to the left 
as e ideeanaes: In practical terms, when the source is closer 
to the receiver than the target is to the receiver, that is, 
when S is less than R, the error contribution of the y value | 
is relatively insignificant. If the geometry is such chat e 
is large, then the error would be paiueed if the angle y 
remains greater than 90°. All plots’ shown in the figures | 


represent the absolute values of the error coefficients. 
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Figure 6 Error vs. Angle for e=0.2 | 
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| Figure 7 Error vs. Angle for e=0.3 
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Figure 8 Error vs. Angle for e=0.4 
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Figure 9 Error vs. Angle for e=0.5 
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Figure 10 Error vs. Angle for e=0.6 
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B. VERIFICATION OF MODEL 
Verification of this model was done by graphical 


confirmation (elliptical) of theoretical expressions. The 


three variables that contribute errors to a particular 


bistatic geometry are t,, t,, and y. The variables t, and ts 


Getermines the eccentricity of the ellipse and y determines 


the jicea tion: of the target on spat ellipse. For each of the. 


three curves, one of the theese variables is assumed to be in. 


error by approximately 10% while the remaining two variables 


were considered to be exactly determined. This forces the two 


remaining terms to become zero, and the fraction«1l error, 
AR/R, is now equal to the one remaining term. The fractional 
error can be crudely determined from the elliptical figure at 
various values of y and compared with the value of ‘the 
appropriate coefficient A,, A,, or Ag ‘obtained from one of the 
figures (Figures 5 - 13) fora corresponding value of e. The 
first of the terns to be validated is the A, term. © Two 
ellipses are drawn in Figure 14 that share the same focal 
points, because t, is held constant, | thus guaranteeing a 
constant inter-focal distance s. Ellipse — 1 has an 


eccentricity of 0.80. If the t value is increased by 10%, 


the eccentricity can be seen from 











ie 














to be reduced to 0.73. Ellipse 2 is the ellipse that peedite. 
from this reduced eccentricity with the same source to 
receiver spacing. At any given value of y, there are two 
corresponding values of R, and the difference between these 
two is AR. This value is divided by the mean of the two R's, 


to give AR/R. In Figure 14 there are four comparisons made at 


separate values of y. There is an illustration of what the 


ae 


RR 





“th een, 


Ellipse 2 
e=0.73 


Ellipse 1 
e=0.80 


Figure 14 Validation of A, 
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a small value of y, the AR/R should be relatively small since 


R has its largest values here. As y increases, AR increases 


as R decreases resulting in a rapid climb of AR/R which is 


‘confirmed in Figures'll and 12. Finally, as y approaches 


180°, AR and R remain constant which explains the asymptotic 
behavior of the function at higher values of y. | 

The next term, A,, is validated through Figure 15, where 
t, and y are held constant and t, is decreased by 10%. Ellipse 


1 in Figure 15 is the original geometry for a given t, and 


J Aa 





_.. ¥igure 15° validation of A, 


ellipse 2 shows the geometry with t, reduced by 10%, which 
gives ab reduced - eccentricity. This graphic is most . 


interesting because there is an immediate explanation for the 
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zero value of A, at about 20° in Figure 11. Shortly beyond 


this value of y, the plot behaves quite similarly to the A, 
plot. As y increases beyond about 90°, the graph again 
demonstrates asymptotic behavior. | 


The validation of the A, function is shown in figure 16. 


AR 
Via 





e=0.80 


Figure 16 Validation of A. 


The validation in this case requires that t, and t, remain 
constant while y is permitted to change slightly. Since t, 
and t. remain constant, so must e remain the same. Thus there 
is only one ellipse to be considered here. Figure 16 shows 
four angles y and four associated Ay's shown by the dotted 
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lines. The difference in length between the dotted line and 


the solid line is the AR value. Note that AR is zero where y+ 


is either O° or 180°. The graphic demonstrates that for 
higher values of e, as in this case, there will be a maximum 
value of AR between 20° and 50°, followed by a gradual deci ine: 
to zero as y increases to 180°. For lower values of 
eccentricity, the AR value reaches its maximum close to 90° 
and has a much smaller value than in the case of larger 
eccentricity. 

The geometrical analysis resulting in Figures 14 through 
16 gives a graphical interpretation of the errors arrived at 
algebraically in Chapter III. From these figures one can 
readily see which errors are more significant for the various 


conditions that may occur in the range problen. 
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V. EXAMPLE CASES 


It is helpful to illustrate use of the analysis by 
applying it to practical problems. The following examples 
will demonstrate the previously mentioned error propagation in 
cases that might be encountered at sea. These examples will 
be approached as if the calculation were conducted from input 
provided by Sonar Control onboard an ASW ship. The 
uncertainty numbers used in these examples are not necessarily 
realistic, but are useful in the demonstration of the method. 
The first three examples are of geometries of low eccentricity 


while the remaining examples are of high eccentricity. 


A. EXAMPLE 1: LOW ECCENTRICITY; NO UNCERTAINTIES IN V 
For the first example assume that the mean sound speed is 
the same over both paths, that is, U equals V. Figure 17 is 
a diagram of an elliptical arrangement of low eccentricity; Ss 
is less than X. Suppose that the following data are provided: 
Y = 110° + 1° 
. tg = 42.62 + 0.018 
Cy = 213.9 + 0.033 
V = 4990 ft/s 
U=V | 
conversion. = 6076 ft/nm 
Calculated values of the various terms in the equations are 


then as follows: 
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C 
Cy 


V(t/-c,? : 
p= —Vltsnts) = 78.966 nm 
2 (t,- t,cosy) 





| ='0.000140 





Ate 
= 0.000235 
C> ! 


Ay = 0.0174 rad 


Target 





42.62 sec Receiver 


Source 
Figure 17 Example Case, Low Eccentricity 
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2 1 
Ay = | — -—————— | = 1.147 
3 i sl 


Ae = 26, _ Cosy _ = -Q.737 
? 1-e2 1-ecosy 


A, = —<&Si2Y = -9.176 


¢  1-ecosy 
If it is assumed that the uncertainties are random 
(independent) then the RMS fractional error may be calculated 
from (3.11) in which it is assumed that U = V, so that 1, = ty, 


and also that AV/V equals 0: 








2 2 
ar| ad ae ad +A? (Ay)? =0.003075 
x x 


Thus the expected uncertainty in range, AR, should be about + 


500 yards or less. 


B. EXAMPLE 2: LOW ECCENTRICITY; UNCERTAINTY IN V 

Assume the same situation as in Example 1, but iow assume 
an uncertainty in Vv of + 3 ft/s. Now, we have an AV/V value 
of 0.0006 in addition to the other errors in Exauple.1. The 


term AV/V in (3.11) is now not zero, but 1, = t, as before. 


‘The fractional error is now — 








| 2 2 - 
(=7y +a ad +A, “e) +A (Ay) 220.,.003133 


Note that the error in AvV/V in this case does not have a large 





(AR) zm8 . 











by le 


effect. 
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C. EXAMPLE 3: LOW ECCENTRICITY; U AND V UNEQUAL 


In this example the situation is the same as in Example 1 
but now let U = 5040 ft/s so that U is not equal to V. The 


data to be used are in this case: 


y = 110° + 1° 
tg = 42.62 + 0.018 
ty, = 213.9 + 0.038 


V = 4990 ft/s 
U = 5040 ft/s 


Now, it is necessary to rescale t, as discussed in Chapter 


II,C., (2.22): 
U ee ge 
ty * Ger = 1.010t, 
From this point, analysis proceeds as in Example 1. 


t 
e= — 2 =0,198 
ty 


V(t,?-t,7) 


oe 2 (t,- &,cosy) 


= 79.873 nm 


2 (2 1 
Rw eo) a 4 
si — ay oe 


“Al. 























As =( -2e ,__ cosy _ = -0.733 


1-e* 1-ecosy 


A; = —7@ SinY 2 _9.174 


1-ecosy 


Now, we use (3.11) with 1, » t,, and AV/V = 0, to get 


2 
=) +A,7 (Ay)? =0.003133 
x  « 





1 | At,\? 
— (AR) .., =,|A (2) +A ( 
?R rms x c. | an 


The mean square fractional error is now somewhat greater than 


that in Example 1, as might be expected. 


D. EXAMPLE 4: HIGH ECCENTRICITY; NO UNCERTAINTIES IN V 
Assume that the mean sound speed is the same ever both 
paths. Figure 18 is a diagram of an elliptical arrangement of 
high eccentricity, which means that the source and receiver 
are generally farther apart than the target and receiver. 
Suppose that the following data are provided: | 


: Y = 110°. + 1° 
t, = 85.23 + 0.018 
t, = 106.54 40.033 
V = 4990 ft/s 
U=V 
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V(t,?-¢,7) 
2 (ty - Cgcosy) 


R= = 12.37 nm 





Source. 85.23 sec Receiver 


Figure 18 Example Case, High Eccentricity 





2 1. 
Ax Fe sia) aie 


A, * 22, _ COSY _| . 4.713 
1-e* 1-ecosy 


A> = ze siny . 9 599 
1l-ecosy 


Equation (3.11) gives for this case: 


eT 


Note that inte is: atiiowt a factor of 3 greater than in the 











An) ) +A 2 (Ay)? #0. 010398 


previous examples of low eccentricity. However, for the 
' geometry, the actual ange gerne is 260 yards, about half of 
the expected error of the earlier cases. It is important to 
point out that the fractional error AR/R is a fairly sensitive 
fuction of eccentricity, but the actual zange error depends on 


the size of the ellipse. 
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S. EXAMPLE 5: UNCERTAINTY IN V 
| Here the same data are provided as in Example 4, but in 
this case V has uncertainty of + 3 ft/s. This time we have 
AV/V = 0.0006 to be included in the RMS calulation, as was 
used in Example 2, and again 1, = t,. The RMS fractional 


uncertainty is thus 








2 At,\* At,\* | | 
The expected range error in this case is also about 260 yards. 


¥. EXAMPLE 6: U AND V UNEQUAL 
In this example, as in Example 3, assume that the same 
data are provided as in the two previous examples except that 
U is not equal to V and there is no uncertainty in Vv. Thus, 
YY = 110° + 1° | 
ts = 85.23 + 0.018 
ty = 106.54 + 0.0338 


V = 4990 ft/s 
U = 5040 ft/s : 


Again we must rescale t, by using 1,: 
; | 
and the process proceeds as in Example 4. 
es ts * 0.792 
ty 


Use (3.11) with t, » 1,, and AV/V = 0 to get a fractional error 
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V(t,2-t,2) © 


——*_S"__ = 12.959 nm 
2 (ty ~ t,cosy) 


R= 





2 1 
A, = -——_+_| = 4, 
- = ae 20 


A. = {222 4-COSY_| = - 4.5290 
d 1-e2 1l-ecosy}/ — 


Ag = ——S2BY = -0.586 
l-ecosy _ 


At, 


2 
1, At 
= (AR) - an( <=) +A { 3 
R = xT .. 3 


2 | 
+A? (Ay)? =0.010327 
Ty oa. 





The expected range error is here about 270 yards, slightly 


larger than in Examples 4 and 5. 
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VI. CONCLUSION 


A. DISCUSSION AND RESULTS 

Improved quieting techniques by submarines have been very 
effective in reducing the exploitable emissions for passive 
sensors. On the other hand, detection by active sonar is 
independent of the submarine’s quieting efforts and gives the 
added bonus of providing an immediate determination of range 
to the target. Thus, active sonar seems to give certain 


obvious advantages. Additionally, a bistatic sonar technique 


increases detection range of active sonar and reduces the E 


vulnerability of the active platform. 

This thesis has examined a method of predicting range and 
range uncertainties in bistatic sonar operations. The method 
described here offers flexibility over conventional methods 
that rely heavily on communication between source and receiver 


for successful operations. Operational methods that use the 


separation distance S are often preferred over those that 


| ‘andice a knowledge of the mean sound speed V, especially near 
fronts and eddies where the mean sound speed may change 
drastically in a few miles. However, this thesis has shown 
that the preference of using the baseline over the mean sound 
speed is unfounded. | _ | 

As shown in the. previous. chapter, the range error is 
controlled to a great extent by the eccentricity of the range 


ellipse, the largest errors occurring in the cases of high 
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eccentricity. Examples are given to provide representative 
applications of the range Equation (2.3) and <5 give the 
estimates of the range errors. | 

We have seen that fractional error AR/R tends to be smail 
when e is small, although dependance an - is still quite 
significant. This suggests that small base-line geometries 


(small values of S leading to small e) are to be preferred if 


' these are operationally feasible. 
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